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Abstract 

We consider a rational system of first order difference equations in the plane with four 
parameters such that all fractions have a common denominator. We study, for the different 
values of the parameters, the global and local properties of the system. In particular, we 
discuss the boundedness and the asymptotic behavior of the solutions, the existence of 
periodic solutions and the stability of equilibria. 
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1 Introduction 

In recent years, rational difference equations have attracted the attention of many researchers 
for varied reasons. On the one hand, they provide examples of non-linear equations which 
are, in some cases, treatable but whose dynamics present some new features with respect to 
the linear case. On the other hand, rational equations frequently appear in some biological 
models and, hence, their study is of interest also due to their applications. A good example of 
both facts are Ricatti difference equations; the richness of the dynamics of Ricatti equations 
is very well-known (see, for instance, [7] and [8]) and a particular case of these equations 
provides the classical Beverton-Holt model on the dynamics of exploited fish populations [3] . 
Obviously, higher order rational difference equations and systems of rational equations have 
also been widely studied but still have many aspects to be investigated. The reader can find 
in the following books [T] . [5] . [TTj . and the works cited therein many results, applications and 
open problems on higher order equations and rational systems. 

A preliminar study of planar rational systems in the large can be found in the paper [B] 
by Camouzis et al. In such work, they give some results and provide some open questions 
for systems of equations of the type 

Q-i + f3ix n + 7iy n 



X ' l+1 Ai + B x x n + C x y n 

Q2 + foXn + 72 Vn 

A 2 + B 2 x n + C 2 y n 



0,1,.. 
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where the parameters are taken to be non-negative. As shown in the cited paper, some of 
those systems can be reduced to some Ricatti equations or to some previously studied second- 
order rational equations. Further, since for some choices of the parameters one obtains a 
system which is equivalent to the case with some other parameters, Camouzis et al. arrived 
at a list of 325 non-equivalent systems to which it should be focused the attention. They list 
such systems as pairs (k, I) where k and I make reference to the number of the corresponding 
equation in their Tables 3 and 4. 

In this paper, we deal with the rational system labelled as (21,23) in [BJ. Note that 
for non-negative coefficients such system is neither cooperative nor competitive but it has 
the particularity that denominators in both equations are equal. This allows us to use 
some of the techniques developed in [2] to completely obtain the solutions and give a nice 
description of the dynamics of the system. In principle, we will not restrict ourselves to the 
case of non-negative parameters, although this case will be considered in detail in the last 
section. Hence, we will study the general case of the system 



where the parameters a\, oli, fix, /?2 are given real numbers and the initial condition (xq, yo) 
is an arbitrary vector of R 2 . It should be noticed that when aifc = otifix the system 
can be reduced to a Ricatti equation (or it does not admit any complete solution, which 
occurs for = f$2 — 0) and therefore these cases will be neglected. Since we will not 
assume non-negativeness for neither the coefficients nor the initial conditions, a forbidden 
set will appear. We will give an explicit characterization of the forbidden set in each case. 
Obviously, all the results concerning solutions that we will state in the paper are to be apply 
only to complete orbits. We will focus our attention in three aspects of the dynamics of the 
system: the boundedness character and asymptotic behavior of its solutions, the existence 
of periodic orbits (and, in particular, of prime period-two solutions) and the stability of the 
equilibrium points. It should be remarked that, depending on the parameters, they may 
appear asymptotically stable fixed points, stable but not asymptotically stable fixed points, 
non-attracting unstable fixed points and attracting unstable fixed points. 

The paper is organized, besides this introduction, in three sections. Section 2 is devoted 
to some preliminaries and some results which can be mainly deduced from the general 
situation studied in [2]. Next, we study the case /?2 = since such assumption yields 
the uncoupled globally 2-periodic equation y n +i — cx-i/Vn and the system is reduced to a 
linear first order equation with 2-periodic coefficients; this will be our section 3 below. The 
main section of the paper is section 4, where we give the solutions to the system and the 
description of the dynamics in the general case ^ 0. We finish the paper by describing 
the dynamics in the particular case where the coefficients and the initial conditions are taken 
to be non-negative. 




71 = 0,1,... 



(1) 
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2 Preliminaries and first results 



Systems of linear fractional difference equations X n+ \ = F(X n ) in which denominators are 
common for all the components of F have been studied in [2] . If one denotes by q the mapping 
given byg(oi,02,...,Ofc+i) = (a 1 /a k+ i,a 2 /a k+ i, . . . ,a k /a k+1 ) for (oi, a 2 , ■ ■ ■ , a k +i) G M fc+1 
with a k+ i ^ and £ : lR fc — > R k+1 is given by £(ai, 0,2, ■ • • , a k ) — {a\, a 2 , ■ ■ ■ , a k , 1), it is 
shown in such work that the system can be written in the form X n+ \ = qo Ao £[X n ), where 
A is a (k + 1) x (k + 1) square matrix constructed with the coefficients of the system. In the 
special case of our system ([!]) one actually has 



Xn+l 
Vn+l 



= qo 




This form of the system let us completely determine its solutions in terms of the powers of 
the associated matrix 

f Pi ai\ 
A = \ h a 2 . (2) 

V 1 J 

Actually, the explicit solution to the system with initial condition (xo,yo) is given by 

(x n+1 ,y n+i y =qoA n (x ,y ,l) t , (3) 

where M* stands for the transposed of a matrix M. Therefore, our system can be completely 
solved and the solution starting at (xq, yo) is just the projection by q of the solution of the 
linear system X n+ \ = AX n with initial condition Xq = (xq, yo, 1)' whenever such projection 
exists. 

Remark 1 When such projection does not exist, then (xo,yo) ^ es ^ n the forbidden set. 
Clearly, this may only happen when for some n > 1 one has 

(0,0, l)A n (x , y , 1)* = 0. 

Therefore, if ai(n) 6 1, < i < 2 are such that A n = ao(n)I + ai(n)A + a 2 (n)A 2 , then one 
obtains immediately that the forbidden set is given by the following union of lines 

F = (J { (x ,yo) G K 2 : ai(n)y Q + a 2 (n)/3 2 x + a 2 (n)a 2 + a (n) = }. 

n>l 

The explicit calculation of ai(n), < i < 2 for each n > 3 may be done in several ways. 
For instance, one has that ao(n) + a\{n)x + a 2 (n)x 2 is the remainder of the division of x n by 
the characteristic polynomial of A. Further, by elementary techniques of Linear Algebra one 
can also compute them in terms of the eigenvalues of A (an approach using the solutions to 
an associated linear difference equation may be seen in 

Remark 2 As mentioned in the introduction, all through the paper we will consider that 

fcctx f Pia 2 . (4) 
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(this is to say, that the matrix A is non-singular) since the cases with foai — /3i«2 may 
be reduced to a single Ricatti equation. Actually, if 012 — P2 — 0, then the system does not 
admit any complete solution, whereas, for 012 7^ or $2 7^ 0, one has that there exists a 
constant C such that ct\ = Cai and (3\ = C/?2 and hence the first equation of the system 
may be substituted by x n +i = Cy n +i and then the second one reduces to the Ricatti equation 

a 2 + /3 2 Cy n+ i 
y n +2 = , n = 0, 1, . . . 

Vn+l 

with initial condition Ui = a2+ P 2X ° , 
a yo 

Our main goal will be to give a description of the dynamics of the system in terms of 
the eigenvalues of the associated matrix A given in ([2]) . We begin with the following result 
concerning 2-periodic solutions which is the particularization to our system of the analogous 
general result given in Theorem 3.1 and Remark 3.1 of [2]. 

Proposition 3 Let us consider the system fip with ai/?2 7^ oli^x- One has: 

1. If fj 2 7^ 0, then there are exactly as many equilibria as distinct real eigenvalues of 
the matrix A. More concretely, for each real eigenvalue A one gets the equilibrium 



( 



2. When f3 2 = 0, one has: 

(a) if ct2 < 0, then there are no fixed points; 

(b) ifO < 0L2 7^ 0\, then there are two fixed points at ( y'ajj o,nd ^-7=^^-, — ^Ja-^j 

(c) if 012 = (i\ and ct\ 7^ 0, then the only equilibrium point is 1 ~ flij / 

(d) if a.2 = Px and ot\ = 0, then there is an isolated fixed point (0, — fix) and a whole 
line of equilibria (xq,(3i). 

3. There exist periodic solutions of prime period 2 if and only if ' ctifa = 0. 

Proof. As stated in [2], a point (a, b) G K 2 is an equilibrium if and only if (a, 6,1) is 
an eigenvector of the associated matrix A. When /?2 7^ 0, it is straightforward to prove 
that for each real eigenvalue A the vector ^ A , A, lj is an eigenvector. In the case 

[3 2 = 0, the equilibrium points can be easily computed directly from the equations 012 = y 2 , 
ai + Pix = xy. 

For the proof of affirmation (3) it suffices to bear in mind that, according to [2], the 
existence of prime period-two solutions is only possible when the associated matrix A has 
an eigenvalue A such that —A is also an eigenvalue. Since A is a 3 x 3 square matrix, this 
obviously implies that the trace of A is also an eigenvalue. Hence, /?i is an eigenvalue but 
this is only possible if ai/32 = 0. If ot\ = 0, then the initial condition (0,yo) gives a prime 
period 2 solution whenever y^ 7^ 122, whereas, if ot\ 7^ and /?2 = 0, a direct calculation 
shows that the solution with initial conditions (0, — /3i) is periodic of prime period 2. □ 

We now study the stability of fixed points in some of the cases. Recall that a fixed point 
of our system (x* , y*) always verifies y* = A for some real eigenvalue A of the matriz A. We 
will say in such case that the fixed point (x*,y*) is associated to A. 
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Proposition 4 Let us consider the system with a.\p2 ^ oli&x- Let p{A) be the spectral 
radius of the matrix A given in 0j and let A be an eigenvalue of A. 

1. If |A| < p(A), then the associated equilibrium is unstable. 

2. If |A| = p(A) and all the eigenvalues of A whose modulus is p(A) are simple, then 
the associated fixed point is stable. Further, if in this case A is the unique eigenvalue 
whose modulus is p(A), then it is assymptotically stable. 

Proof. The Jacobian matrix of the map F(x,y) = ( y ai+ y lX , a2+ f 2X ^j at a fixed point 
(x* , y*) is given by 

' fc/y* -x*/y* ' 



DF(x*,y*) 



Nv* -i 



Consider an eigenvalue A of A and let A2, A3 be the other (non- necessarily different) eigenval- 
ues of A. Let us show that the eigenvalues of the Jacobian matrix at a fixed point associated 
to A are just ^ and The result is trivial when /?2 = since the eigenvalues of A are ft\ 
and ± v /«2 an d fixed points are always associated to one of the eigenvalues ±^/o~2. If P2 7^ 0, 
then x* — — -p^ 2 - and y* = A and, therefore, one obtains: 

pi — A A2 + A3 



trace(DF(x*,y*)) 
det(DF(x*,y*)) 



A A 

~/3iA + A 2 -a 2 det(A) A 2 A 3 



A 2 A 3 A 2 



showing that the eigenvalues of DF(x* ,y*) are as claimed. Now, the first statement follows 
at once since, if |A| < p(A), then at least one of the eigenvalues of DF(x* ,y*) lies outside the 
unit circle. Moreover, when |A| = p{A) and it is the unique eigenvalue with such property, 
then the eigenvalues of DF(x* ,y*) are inside the (open) unit ball and hence the equilibrium 
(x*,y*) is assymptotically stable, which proves the second part of (2). 

For the proof of the first part of (2) let us recall that if (x*,y*) is a fixed point of ([T|) 
associated to the real eigenvalue A, then X* = (x* , y* , 1)' is a fixed point of the linear system 
X n +i — jAX n . The eigenvalues of the matrix M — jA are obviously 1, ^ and Since 
the eigenvalues of A having modulus p(A) are simple, so are the eigenvalues of M having 
modulus 1. Therefore, the fixed point X* is stable [51 Th. 4.13]. Now, the stability of 
(x*,y*) follows at once from and the continuity of q in the semi-space z > 0. □ 



3 Case (3 2 = 0. 

Recall that, since we are assuming that inequality (U]) holds, we have f3ia.2 ^ 0. In this case, 
the forbidden set of the system reduces to the line y = 0. Since f3 2 — the second equation 
of the system becomes the uncoupled equation 

Ci2 

Vn+l = , 

Vn 
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which, as far as a,i ^ 0, for each initial condition yo ^ gives 

yo for even n, 

0.2 

— for odd n. 

Vo 



(5) 



Substituting such values in the first equation of the system we obtain a first order linear 
difference equation with 2-periodic coefficients whose solution is given by x\ = (cti+ 0xo) / yo 
and, for n > 1, 



a 2 



x 



ai(/3i + y ) 



Ml | ( 

yo yo V "2 



a 2 



a 2 



k=l vpi 



Q'2 



fc=l VPl 



for even n, 



for odd 



(0) 



Hence, we have proved the following. 



Proposition 5 If /?2 = and /?2«i ^ /?ia2; i/ien £/ie system (QJ) is solvable for any initial 
condition (xo,yo) with yo ^ and t/ie solution (x n ,y n ) is given by |5[) and (0|) where, ex- 
plicitly, one has: 



(1) If a 2 = 01, then for n > 1 



< 



x 



+ yo ) n 
2/3i 2 



ai A .tq _ aiQfli + y )(re - 1) 
L yo yo 2/3iy 



/or even re, 
/or odd re. 



(2) Ifa 2 + 0{, then for re > 1 



_L 

a- o 



2\ 2 



X 



ai(/3i + y ) 
fii - a 2 



1 - 



C>2 



yo yo V "2 



Xo 



ai(0 + yo) 



i - 



OL2 



a 2 

A 3 



for even n, 



for odd re. 



From the proposition above one can easily derive the following result which completely 
describes the asymptotic behaviour of the solutions to the system. 

Corollary 6 Let us consider 02 — and 0ct2 7^ 0. 

1. When fti 2 — a.2 one has: 

(a) Ifaij^O, then every solution to the system is unbounded except those with initial 
condition (xq,—0i), which are 2-periodic. 
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(b) If a± = 0, the system is globally 2-periodic. 

2. If Pi 2 = —ct 2 , then the system {IP is globally J^-periodic. Further, the solution corr 
sponding with the initial condition (xo,J/o) * s of prime period 2 if and only if2p 2 x Q 
aiWi+yo) =0. 



3. If Pi 2 7^ I c^2 1 7 then the solutions with initial condition [ — — — — —Sr-^- , yo I are period- 

V a 2 -Pi J 

two solutions. Moreover, 

(a) If (3 1 2 > [as], then any other solution to the system (QP is unbounded. 

(b) If fti 2 < |ct2 1 , then any other solution of $7$) is bounded and tends to one of the 
period-two solutions described above. 

Proof. The proof is a straightforward consequence of the explicit formulas for x n and y n 
given in Proposition [SJ It should, however, be mentioned that the globally periodicity of the 
system in the case p\ 2 = —a-i can be easily seen since the associated matrix A given by @ 
in such case verifies A 4 — f3fl, where I stands for the identity matrix. Actually, a simple 
calculation proves that the solution starting at (xo,yo) is the 4-cycle 

(ai + PxXQ -f3 2 \ ( ai(pi+y Q ) \ f-P 2 x + a ± y ~P\ 
{xo,yo),\ , , \-xo m ,y 



oci (Pi + y ) 



yo yo J \ Pi J \ Pm 2/0 

which is obviously 2-periodic if and only if xq — —xq — ai it .□ 

From the above result and Proposition |4] one easily gets the following information about 
the stability of the fixed points. 

Corollary 7 Let us consider p2 = and P1Q.2 =/= 0. 

1 . If Pi 2 = ct2 , then 

(a) for cti ^ the unique fixed point of {!]) is unstable; 

(b) for cti = every fixed point of £2) is stable but not asymptotically stable. 

2. If Pi ^ a 2 > 0, then 

(a) for Pi 2 > q;2 both fixed points of f7|] are unstable; 

(b) for Pi 2 < ct2 the fixed points of fljj are stable but not asymptotically stable. 

4 Case f3 2 ^ 0. 

Proposition 8 Suppose p2 ^ and (xo,yo) is an initial condition not belonging to the 
forbidden set F. In such case the solution of system {!]) is given by 

_ Vn+l 1 02 _ Vn 

V n -i P 2 P2 Vn-l 

where v n is the unique solution of the linear difference equation 

Vn+3 - PlVn+2 ~ Ci 2 V n+ i + [PlOt 2 ~ p\ot\) V n = (8) 

with initial conditions u_i = 1, vo = yo and vi = p 2 xo + ct 2 . 
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Proof. As we have seen in section 2, the solution to system ([T]) starting at a point {xo,yo) 
not belonging to the forbidden set is just the projection by q of the solution of the linear 
system (ti n +i, tVt+i, ifn+i)* = A(u n ,v n ,w n Y with initial condition (xo,yo,lY where A is 
given by @. Since the third equation of such linear systems reads w n +\ = v n , it can be 
reduced to the planar linear system of second order equations 

Un+l = ft^n + «li;„_l, 

and hence, if (u n , v n ) is the solution to ^ obtained for the initial conditions (uo, vq, ^-i) = 
(a:o, 2/o j 1): then the solution of our rational system for the initial values (xq, yo) will be 

x n+1 = u n /v n -i , y„+i = v n /vn-i- (10) 

It is clear that for /3a ^ we have that u n can be completely determined by ([9]) in terms of 
v n+ \ and u n _i and hence it suffices to solve the third order linear equation 

v n +3 - PiV n +2 - a2V n +i + (Pic® - fhoti) v n = (11) 

trivially deduced from © and substitute the corresponding values in (jTCTj) to obtain the 
result claimed. □ 

In the following results we shall discuss the behavior of the solutions to ([I]) by using 
Proposition [8] We shall consider three different cases depending on the roots of the charac- 
teristic polinomial of the linear equation (j8]). Recall that such roots are also the (possibly 
complex) eigenvalues of the matrix A given in (|2]l. 

From Proposition [5] we see that the asymptotic behavior of the solutions of system (JXJ) 
will depend on the asymptotic behavior of the sequences v Vn , being v n solutions of the 
linear difference equation ©. The Theorem of Poincare (Theorem 8.9 in [5]) establishes a 
general result for the existence of lim^oo rp 24 -. In our case, since equation flSJ) has constant 
coefficients, we can directly do the calculations, even in the cases not covered by the Theorem 
of Poincare, to describe the dynamics of system ([T]). 



4.1 The characteristic polinomial has no distinct roots with the 
same module 

Let Ai, A2, A3 be the three roots of the characteristic polynomial of the linear difference 
equation (jHJ in this case. A condition on the coefficients for this case can be given by: 



with «i ^ or a2 < 0. Recall that we assume here that fiia.\ ^ and (3-2 ^ 0. 
If Ai is the characteristic root of maximal modulus, we will denote by L the line 

L = {(x,y) : fox=(ft 1 -\ 1 )(y + \ 1 ) }. 
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Proposition 9 Suppose that f3 2 ^ and every root of the characteristic polinomial of the 
linear difference equation ftB\) is real and no two distinct roots have the same module. When 
(xo,yo) is not in the forbidden set, we have: 

1. J/|Ai| > |A 2 | > |A 3 |, then 

( A 2 — Ct 2 \ 

(a) the system 0] admits exactly the three equilibria I - , Ai J, i = 1,2,3; 







2 



Ai — «2 



(b) the fixed point ( , Ai ) attracts every complete solution starting on a point 

P2 



(xo,yo) which does not belong to the line L; 

initial rnnd.itinn. ( Tc\ =L I — — 

P2 



' \ 2 — Q, 2 

(c) the corresponding solution to the system with initial condition (xq, yo) 7^ ( t. , A3 



(A 2 ^ — Q?2 
a ' ^ 2 
P2 

2. If I Ai| > J A2 j md Ai has algebraic multiplicity 2, then 



\ 2 

Aj - «2 



(a) the system 0) admits exactly the two equilibria 1 , Ai ], i= 1,2; 

V P2 / 

( Ai — GL2 A 

(b) the fixed point , Ai I attracts every complete solution except the other 

V P2 / 

fixed point; 

3. If I Ai| > |A 2 | and X 2 has algebraic multiplicity 2, then 

(a) the system 0] admits exactly the two equilibria ( 1 — , Ai J , i — 1,2; 



2 



\ 2 

Ai — a-i 



(b) the fixed point , Ai I attracts every complete solution starting on a point 

\ P2 / 
(xo,yo) which does not belong to the line L; 

(c) the corresponding solution to the system with initial condition (xo,yo) G L con- 

( A 2 - «2 , 

verges to , A 2 

V P2 

4- If Ai has multiplicity 3, then 

f\ 1 2 -a 2 \ 

(a) the system flp has a unique equilibrium , Ai I; 

V P2 / 

(b) the equilibrium is a global attractor. 

Proof. In all the cases, the equilibrium points are directly given by Proposition [3] The 
assertions concerning the asymptotic behaviour can be derived as a consequence of Case 1 
in [5J pg. 240], bearing in mind that 

v n +i 1 a 2 v n 

S- > Vn 



U n _l [3 2 P2 V n -l 
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and that v n is the solution to the linear equation (JSJ with initial conditions w_i = 1, vq = ya 
and vi — fox® + a^-Q 

4.2 The characteristic polinomial has two distinct real roots with 
the same module 

It is easy to check that this case occurs when p\ ^ 0, (3 2 ^ 0, a\ = and a2 > 0. Thus, the 
roots of the characteristic polynomial of the linear difference equation ([5]) are /3i and ±Ja2~. 



Proposition 10 Suppose (3\ ^ 0, /?2 ^ 0, ct\ — and o?2 > 0. Assume also that (xo,yo) is 
not in the forbidden set. 

1. If Pi 2 — a 2 , then 

(a) there are two equilibrium points (0,±/3i); 

(b) the equilibrium point (0, attracts every complete solution not starting on a 
point of the line x = 0; 

(c) the solutions starting on a point (xo,yo) of the line x = are prime period-two 
solutions except the two equilibrium points (0,±p\). 

2. If Pi 2 > a 2 , then 

fPi 2 -a 2 \ 

(a) there are three equilibrium points , (3\ I and (0, i^/a^); 



Ih 



a 2 

Pi - a 2 



(b ) the equilibrium point I , fi\ I attracts every complete solution not starting 

V P2 / 

on a point of the line x = 0; 

(c) the solutions starting on a point (xo,yo) of the line x = are prime period-two 
solutions except the two equilibrium points (0, ±y / a2); 

3. If (3 2 <a 2 , then 



(p 2 -a 2 \ 

(a) there are three equilibrium points ,/3i I and (0, iy/a^); 

V P2 / 

(b) the solutions starting on a point of the line x — are prime period-two solutions 
except the two equilibrium points (0, i^/aj); 

(c) the solutions starting on a point of the lines p^H — 2 J^ 1 y = or x = 1 - 

Pi ' P2 

are unbounded with the only exception of the fixed point ( , j3\ ] ; 

V P2 / 

(d) the solutions starting on any other point (xo,yo) are bounded and each tends to 
one of the two-periodic solutions; 
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Proof. In all cases the affirmation a) is a consequence of Proposition [3] 



When Pi 2 = a 2 , the roots are /3i, with algebraic multiplicity two, and —f)\. By Propo- 
sition [5] we know that any solution of the system can be written as 

_ (n+l)P 1+ P 2 + P 3 (-l)" +1 2 R 2 
nPi+P 2 + P 3 (-l)" 

where Pi, P 2 and P3 actually satisfy 

Pi+P 2 -P 3 = /32X ° +/?1 " , P 2 + P 3 =2/o, -Pi+P 2 -P 3 = /3i. (12) 
Pi 

If Pi 7^ 0, then (x n ,y n ) obviously tends to (0,/3i). From (fT2|) we see that Pi = if 
and only if xo — and, in such = and y n takes alternatively the values Afii and 

A~ x f3i with A = p^p^ . Notice that t/o^0 guaranties P 2 + P 3 7^ and, since p\ 7^ 0, we 
can not have Pi — and P 2 — P3 = 0. This completes the proof of (1). 

In the case f3f 7^ a 2 , by Proposition [8] we can write the general solution of the system as 



, n+l 

2 



P 1+ [P 2 + P 3 (-1)" +1 ](^) 

Pl%n = T^=TPl - "2, 

Pl + lP.+P^-l)- 1 ]^, 

Pl + [P 2 +P 3 (-l) r 
Vn = 



V a 2 



n-1 • 



Pi + IP^ + p^-i)"- 1 ] (^r 

where Pi, P 2 and P3 satisfy 

PiA + (P 2 - P 3 )V^2 = /3 2 so + a 2 , 

Pi+P 2 +P 3 = 2/o, (13) 

Pi/3r 1 + (p 2 -p 3 )V^ T - 1. 

When /3f > a 2 one immediately gets the results of statement (2) with an argument similar 
to that of the previous case. Therefore we will focus our attention in the case /3f < a.%. The 
condition x — is, according to (fT3|) . equivalent to Pi = and in such case one gets 
x n — and y n takes alternatively the values K^/cn and K~ x ^faJi with K = pHrp^ = 
Now, if Pi ^ and the initial conditions are taken such that P 2 + P3 7^ 7^ P 2 — P3, then 
(xn,y n ) tends obviously to the 2-cycle {(0, K^/a^), (0, K" 1 ^/^)} where K = pji*p ■ On 
the contrary, if either P 2 + P3 = or P 2 — P3 = (and only one of both equalities holds) then 

both sequences x n and y n are unbounded. From the system (|13|) one gets that P 2 — P3 = if 

a 2 

and only if xq = {01 — ct 2 )//3 2 and that P 2 + P3 = is equivalent to /3 2 xo H 3 j/o = 0. 

Pi 

This shows the validity of c).0 
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4.3 The characteristic polinomial has complex roots 

Now we consider the case in which the characteristic polynomial of the linear difference 
equation has a couple of complex roots pe , with sin 9 > 0. Let A ^ be the real root. It 
can be easily shown that 

(3\ = \ -\- 2p cos 9; a 2 = - (2Apcos6> + p 2 ) ; faai = Ap 2 + /?ia 2 . (14) 

and that this situation occurs when 



|gl02 - ftai - ^/?i 3 ' 



2 




(\ 2 ~a 2 \ 

By Proposition [3] we know that the unique equilibrium is , A I . Denote by L 

V P2 / 

the line 

L = { (x,y) : /3 2 x = (/3 1 -\)(y + \) }. 

Notice that (/3i — A)(y + A) = 2ypcos9 — a 2 — p 2 . Also observe that the equilibrium does 
not belong to L. 

Theorem 11 Suppose f3 2 ^ and the characteristic polynomial of the linear difference 
equation to have complex roots and assume that (xo,yo) is not in the forbidden set. 

1. The solutions starting on the line L remain on it and they are either all periodic or 
all unbounded. 

2. If |A| > p, then the unique equilibrium attracts all the solutions not starting on L. 

3. If |A| < p, then every non fixed bounded subsequence of a solution accumulates on L. 

4- If |A| = p, then every complete solution (neither starting on the fixed point nor on L) 
lies on a not degenerate conic, which does not contain the equilibrium. 

Proof. Assume that (£0,2/0) is not the fixed point. Using Proposition El we have 

PA" +1 + 2p n+1 cos (a+(n+l) 9) 



a 2 + (3 2 x n = 



y„ 



PA"" 1 + 2p n ~ 1 cos (a + (n- 1) 9) 
PA" + 2p n cos (a + n9) 



PA"- 1 + 2/?"- 1 cos (a + (n - 1) 9) ' 
where the constants P £ K and a 6 [0, 2ir), together with k £ M + , are given by 
A pe lB pe- 4 ' 

1 1 111 fee- 1 = 1 y I . (15) 





1/A e- w /p e w /p 

Observe that we may consider P > 0, by replacing, if necessary, a with a + it. 
Let us consider the sequences 

a n = 2 y— J cos (a + nO) ; r„ = 2 (— J sin (a + n9) . 
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It can be easily proved that 

2 P + (T n+1 P + (T n 

a2 + P2x n = X — , y n = A—— , (16) 

Acr n +i = pa n cos — pr n sin 0, pc n -i = Xcr n cos + Ar^ sin . (17) 

As a consequence, A 2 ct„+i — 2Xpa n cos0 + p 2 a n -i — 0, and then 

2 „A 2 -2pAcos6> + /9 2 

a 2 + Pix n = 2py n cos - p + P — , 

P + a n -i 

which is equivalent to 

„A 2 — 2pAcos6> + p 2 , s 

fox n -(f3 1 -X)(y n + X)=P -p- F -. (18) 

" + 0"n-l 

Using (|15p one has that (xo,yo) € L if and only if P = and from (|18[) we then get that 
{x n , Vn) G L for all n > 1. 

Furthermore, by (fTB"]) . we see that if (xo,yo) G A then the solution (x n ,y n ) is periodic 
whenever ^/V is a rational number and unbounded otherwise. 

Assume now that the solution (x n ,y n ) does not start on L, this to say, P ^ 0. We will 
now distinguish the three cases: |A| > p, |A| < p and |A| = p. 

If |A| > p, then by §TE\i one immediately has x n — > A and y„ — > A. 

Suppose now that | A| < p. If (x„ fc , y„ fc ) is a subsequence satisfying that inf |cos (a + (rik — 1) 0) \ > 

k 

0, then one obviously has <J nk -\ — > oo. Using the definition of a n , one easily gets that 



is bounded. Then, {x nk ,y nk ) is a bounded subsequence and equation (|18j) shows that it is 
attracted by the line L. 

On the other hand, if cos (a + (rife — 1) 0) — > 0, then the left equation in (|17p lead us to 

— ) — > 2 sin# > 0. Thus, er„ fc — » oo and using (TITf once more we get — — > oo. 

Therefore, (x nk ,y nk ) is an unbounded subsequence. 

Finally, let us suppose p = |A|. If we consider the change of variables 



2 \ A , w pXcos0 



x = (fox + a 2 - p 2 ) — — (y-X) 

v ' 2A cos — 2p X cos — p 

to , 2\ 1 / ,(A + pcos6») 



2 sin sin 

then one may deduce from (jTTJ) that 3^" = — " 1 , y^ = — — - 1 . Therefore, one 

P + <7n-l P + (T n -i 

immediately gets that 

x^ + yf= ^ n {pX? r, , (x^- P X) 2 = P 2 {pX? 



(P + a n ^) 2 ' v " r ' {P + Vn-lY ' 

which clearly shows that (x~^,y^) lies in the conic x 2 + y 2 = -pz(x — pX) 2 , having its focus in 
(0,0), its directrix in the line x = pX and eccentricity 2/ P. Further, one immediately sees 
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that the fixed point — , A I is transformed by the change of variables above in (0, 0) 

\ ft / 
and, hence, it does not belong to the conic. □ 

Remark 12 In the case |A| < p of this last theorem, one might conjecture that every subse- 
quence of a solution (even a non-bounded one) actually approaches the line L; but this is not 
the case. Let us take, for example, the system whith ct\ = = 3,a2 = — 4,/?2 = —10, in 
which the characteristic roots of the associated polynomial are given by A = 1 and \/2e I7r / 4 
and consider the solution starting on (xo,yo) = ( — H/20, 3/2). We then have that a = 0, 
P = 1 and <72+4fc = for all k > 0. One may use equation U8\) to show that all the 
points of the form (^3+4^, 2/3+4/c) lay on the line lOx + 2y + 3 = while the line L is 
given by lOx + 2y + 2 = 0. Note, however, that the subsequences (x4fc,y4fc), (xi+4fc, 2/1+4/c) 
and (a;2+4fe) 2/2+4fe) are att bounded and converge respectively to (—3/5, 2), (—2/5, 1) and 
(—1/5,0), which do belong to L. 

It should also be noticed that the fixed point lays on the line 10x + 2y + 3 = 0. This is 
also the case in the general setting. It follows from H8\) that whenever cr nfc _i = then the 
point (x n -k,y n -k) is on the line containing the fixed point which is parallel to L. 

Remark 13 Notice that, according to the results in when |A| = p and the argument 9 
of the complex root is a rational multiple of tt, the system is globally periodic. 

4.4 Stability of fixed points 

We finish this section with the complete study of the stability of the fixed points in the case 
ft #0. 

Theorem 14 Suppose that f3% 7^ 0, let X be a real eigenvalue of the matrix A given in 
Let ( A ~ 2 " 2 , A^ be the associated fixed point and denote by p(A) the spectral radius of A. 

1. If \\\ < p{A), then the associated fixed point is unstable. 

2. If |A| = p{A), then the associated equilibrium is stable if and only if every eigenvalue 
whose modulus is p(A) is a simple eigenvalue. Moreover, the stability is asymptotic if 
and only if X is a simple eigenvalue and it is the unique eigenvalue of A whose modulus 
is p(A). 

Proof. The first statement was already proved in Proposition |4] Besides, in such Propo- 
sition, we have shown that if every eigenvalue whose modulus is p(A) is simple then the 
associated equilibrium is stable. Let us prove the converse. 

According to the results of the previous subsections, the only cases in which one has a 
non-simple eigenvalue of maximal modulus are the cases treated in Proposition [9] (1) and 
(4) and the first case of Proposition fTOT We will see that in such cases the equilibrium points 
associated to eigenvalues of maximal modulus are unstable. 

We begin with the case of an eigenvalue Ai of maximal modulus with multiplicity 2. 
For each N € N, N > 1 one may consider the solution with initial conditions (xo,yo) — 

( Xl p" 2 — (jy5 ji)ff 2 ; ^1 — jy^+i )- The solution of © in such case is given by v n = X" +1 (N 2 + 
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1 — Nn — N)/(N 2 + 1), which cannot vanish since N > 1. For this solution one has 

A 2 —a 

\vn ~ = | Ai | -/V, proving that the equilibrium ( 1 „ 2 , Ax ) is unstable. 

Similarly, if A has a unique eigenvalue A of multiplicity 3 then, for each N £ N, N ^ 
let us consider (x , yo) = ( A — 2 ^ N2 , A). The corresponding solution to is given by 

u„ = ( N ~ n ~^l — • It is not difficult to see that v n ^ for all n > 1 and then the solution 

to our system (JJ) is complete. Further, since y n = v n /v n -i one gets that \y^ — A| = 2|A|. 
Therefore, the fixed point ( A , A) is not stable. 

When ol\ — 0, /3f = a.^ ^ 0, there are two equilibrium points associated to eigenvalues 
of maximal modulus: (0, ±/3i) . The fixed point (0, ~/3i) is, according to the result of 
Proposition [TU1 unstable since the other equilibrium attracts all the solutions not starting 
on the line x = 0. To see that (0,/3i) is also unstable, let us choose, for each odd N € N, 
the solution starting at (xo,yo) = ( Npl > ffi ) • Then, using equation (|12p and the expression 
for y n given just above such equation, we have v n = /3™ +1 + Pin/3™ if n is even and v n — 
pn+i _|_ _|_ if n i s odg^ where Pi = —fli/N. Since iV is odd, we see that y n exists 
for all n e N and, further, we get that |yjy — = 2|/?i|, which clearly implies that (0,/3i) 
cannot be stable. 

Finally, it only remains to prove that when A has distinct simple eigenvalues whose 
modulus equal p(A) then the fixed point is not asymptotically stable. But this situation can 
only happen if either one has the situation described in Proposition [TT] (4) or the one given 
in Proposition [10] (3). In the case of complex eigenvalues we had seen that all the orbits lie 
on conies not going through the fixed point and, hence, it cannot be asymptotically stable. 
In the other case, it is clear that the fixed points (0,± v / O2) are not attracting since every 
solution starting on the line x = is 2-periodic. □ 

Remark 15 It is interesting to notice that in the three cases in which there is an eigenvalue 
of maximal modulus with multiplicity larger than 1, the corresponding fixed point is attracting 
but unstable. 

5 Non-negative solutions to the system with non-negative 
coefficients 

When the coefficients of our system ((TJ) are non-negative and we restrict ourselves to non- 
negative initial conditions, many of the cases studied in the previous sections cannot appear. 
Further, in such case one may describe which kind of orbits appear and their assymptotic 
behaviour without the previous calculation of the characteristic roots. 

It shoud be noticed that whenever the coefficients in system ([1]) are non-negative and 
«i/?2 7^ 012P1, every initial condition (xo,yo) with xq > 0, yo > gives rise to a complete 
orbit except for 02 = where the condition xq > is also necessary. 

It will be convenient to study independently the case ai/?2 = 0. Next result is a simple 
summary of the results in Section 3 and Proposition 1101 and, hence, we omit its proof. 

Corollary 16 Let us consider that the coefficients in system {Ip are non-negative and 
d\t3i = 7^ a 2 /?i. 
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1. If 02 = 0, one has: 

(a) When a-i < P\ there are no non-negative periodic orbits and all non-negative 
solutions are unbounded, with the only exception of the case u<x = P\, a.\ = 0, 
which is globally 2-periodic. 

(b) When ct2 > P\ there exists a non attractive fixed point ^ ^c^-pi ' anc ^ ^ e 
whole line (012 — Pi)xq = + yo) of 2-periodic solutions. Every other non- 
negative solution is bounded and converges to one of the 2-cycles. 

2. If fii 7^ = ai, then every non-negative solution is bounded and the ones starting in 
the line xq = are 2-periodic. Morevoer, 

(a) When 0,2 < Pi there are two non-negative fixed points: ^ 1 ^" 2 , Pij , which at- 
tracts all non-periodic non-negative solutions, and (0,^/02). 

(b) When 0:2 — Pi there is a unique non-negative equilibrium (0,/3i) which attracts 
all non-periodic non-negative solutions. 

(c) When a 2 > Pi the unique non-negative equilibrium is (0,^/02) which is not an 
attractor. Every non-negative solution converges to one of the periodic solutions. 

The remaining cases are jointly treated in the following result. All the definitions and 
results on non-negative matrices which are used in its proof may be found in [TOl Ch. 8] 

Proposition 17 Suppose that system (QJ) has non-negative coefficients and that otifa 7^ 0. 

1. If 012 ^ or Pi 7^ then there is a unique non-negative (actually, positive) stable 
equilibrium which attracts all non-negative solutions. 

2. If ct2 =01 = 0, the system is globally 3-periodic with a unique equilibrium. 

Proof. Let us consider A as in ^ . A simple calculation shows that (^4 + I) 2 is positive 
and, therefore, A is irreducible. Then the spectral radius p(A) is a strictly positive simple 
eigenvalue of A. 

If there exists another eigenvalue A such that |A| = p(A) then, since A is non-negative 
and irreducible, the eigenvalues of A should be Afe + i = p{A)e lk ' K ^ where k = 0,1,2 and, 
consequently, A 3 = p(A) 3 I. The direct computation of A 3 shows that this is possible if 
and only if CK2 = Pi = and, hence, in that case, the system is 3-periodic and the only 
equilibrium is the one associated to the real eigenvalue p{A). 

In the remaining cases, Ai = p(A) is a dominant eigenvalue and, according to our results 
of Proposition 01 Proposition |H] and Theorem [TT1 the corresponding fixed point is stable and 
attracts all complete solutions except those starting on the line 

L = { (x,y) : j 8 2 a; = G8i-A 1 )(i/ + Ai) }. 

Since Ai is the largest eigenvalue of A, one has that det(A— p,I) < for all p, > Ai. However, 
det(A — Pil) = (X1P2 > 0, showing that Pi < Ai. Thus, for every x > and yo > one 
obtains P2X0 > and (Pi — Xi)(yo + Ai) < 0, which proves that (x , yo) $ L. 

The equilibrium associated to the eigenvalue Ai = p(A) is ( Al ^ 2 ° 2 , Ai^ , which is positive 
since, as before, one sees that det(A — y/a2~I) — aip2 > and hence Ai > y/o2- □ 
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